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Abstract. — This paper is a contribution to the general program of embedding theories of 
dynamical systems. Following our previous work on the Stochastic embedding theory devel- 
oped with S. Darses, we define the fractional embedding of differential operators and ordinary 
differential equations. We construct an operator combining in a symmetric way the left and 
right (Riemann-Liouville) fractional derivatives. For Lagrangian systems, our method pro- 
vide a fractional Euler-Lagrange equation. We prove, developing the corresponding fractional 
calculus of variations, that such equation can be derived via a fractional least-action prin- 
ciple. We then obtain naturally a fractional Noether theorem and a fractional Hamiltonian 
formulation of fractional Lagrangian systems. All these constructions are coherents, i.e. that 
the embedding procedure is compatible with the fractional calculus of variations. We then 
extend our results to cover the Ostrogradski formalism. Using the fractional embedding and 
following a previous work of F. Riewe, we obtain a fractional Ostrogradski formalism which 
allows us to derive non-conservative dynamical systems via a fractional generalized least- 
action principle. We also discuss the Whittaker equation and obtain a fractional Lagrangian 
formulation. Last, we discuss the fractional embedding of continuous Lagrangian systems. 
In particular, we obtain a fractional Lagrangian formulation of the classical fractional wave 
equation introduced by Schneider and Wyss as well as the fractional diffusion equation. 
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"II y a de l'apparence qu'on tirera un jour 
des consequences bien utiles de ces paradoxes^ 1 ) , 
car il n'y a gueres de paradoxes sans utilitc." ( 2 ) 

Leibniz, Letter to L'Hospital, September 30, 1695 

Comme la construction du monde est la plus parfaite possible 
et qu'elle est due a un Createur infiniment sage, 
il n'arrive rien dans le monde qui ne presente quelqucs proprictc dc maximum ou minimum. 

Aussi ne peut-on douter qu'il soit possible de determiner tous les eflets de l'Univers 
par leurs causes finales, a l'aide de la mcthodc des maxima et minima 
avec tout autant de succes que par leurs causes efficientes. 

Methodus inveniendi lineas curveas maximi minimive proprietate gaudentes 

in Eulcrii Opera omnia, Serie I, 24, 
Bcrlin-Bascl-Boston-Stuttgart, Lipsiae-Birkhauscr Vcrlag, 1911 



INTRODUCTION 

The aim of this paper is to introduce a general procedure called the fractional embedding 
procedure, which roughly speaking allows us to associate a fractional analogue of a given 
ordinary differential equation in a more or less canonical way. The fractional embedding 
procedure is part of a global point of view on dynamical systems called embedding theories 
of dynamical systems [15]. Before describing the general strategy underlying all embedding 
theories and the particular fractional embedding procedure, we provide a set of problems 
which lead us to our point of view: 

- Turbulence: Fluid dynamics is modelled by partial differential equations. Solutions 
of these equations must be sufficiently smooth. However, there exits turbulent behavior 
which correspond to very irregular trajectories. If the underlying equation has a physical 
meaning, then one must give a sense to this equation on irregular functions. This remark 
is the starting point of Jean Leray's work on fluid mechanics [29]. He introduces what 
he calls quasi- derivation and the notion of weak-solutions for PDE. This first work has a 
long history and descendence going trough the definition of Laurent Schwartz's distribu- 
tion and Sobolev spaces. We refer to [2] for an overview of Jean Leray's work in this domain. 



W Derivatives of non-integer order 

( 2 )±t will lead to a paradox, from which one day useful consequences will be drawn. 
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- Deformation Quantization problems: The problem here is to go from classical mechan- 
ics to quantum mechanics trough a deformation involving the Planck constant. Roughly 
speaking, we have a one parameter h family of spaces and operators such that they 
reduce to usual spaces and operators when h goes to zero. For example, we can look for 
a deformation of the classical derivative using its algebraic characterization trough the 
Leibniz rule. Another way is following L. Nottale [35] to assume that the space-time at 
the atomic scale is a non-differentiable manifold. In that case, we obtain a one parameter 
smooth deformation of space-time by smoothing at different scales. The main problem is 
then to look for the deformation of the classical derivative during this process. We refer 
to [16] for more details. 

- Long term behavior of the Solar-system: The dynamics of the Solar system is usually 
modelled by a n-body problem. However, the study of the long-term behavior must 
include several perturbations terms, like tidal effects, perturbations due to the oblatness 
of the sun, general relativity effects etc. We do not know the whole set of perturbations 
which can be of importance for the long term dynamics. In particular, it is not clear that 
the remaining perturbations can be modelled using ordinary differential equations. Most 
of stability results uses in the Solar systems dynamics make this assumption implicitly 
[30]. An idea is to try to look for the dynamics of the initial equation on more general 
objects like stochastic processes, by extending the ordinary derivative. Then one can look 
for the stochastic perturbation of the underlying stochastic equation which contains the 
original one. As a consequence, we can provide a set of dynamical behaviors which have 
a strong significance being stable under very general perturbation terms. This strategy is 
developed in ([18] [20]) and applied in [21]. 

These problems although completely distinct have a common core: we need to extend 
the classical derivative to a more general functional space. This extension being given, we 
have a natural, but not always canonical, associated equation. 

However, we are lead to two completely distinct theories depending on the nature of the 
extension we make for the classical derivative. 

For the Solar-system problem, we need to extend the classical derivative to stochastic 
processes by imposing that the new operator reduces to the classical derivative on differen- 
tiable deterministic processes. In that case, the initial equation is present in the extended 
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one and we use the terminology of embedding theory. We describe the strategy with more 
details in the next section. 

For quantization problems, in particular first quantization of classical mechanics, we 
have an extra parameter h (the Planck constant) and the extended operator denoted by 
T>h reduces to the classical derivative when h = 0. The initial equation is not contain in 
the extended one, but we have a continuous deformation of this equation depending on h. 
We then use the terminology of deformation theory. We give more details on this type of 
theories in the following. 



1. Embedding theories 

The general scheme underlying embedding theories of ordinary or partial differential 
equations is the following: 

— Fix a functional space T and a mapping i : C° — > T . 

— Extend the ordinary derivative on T by imposing a gluing rule. 

— Extend differential operators. 

— Extend ordinary or partial differential equations. 

Let us denote by V the extended derivative on T . The gluing rule impose that the extended 
derivative reduces to the ordinary derivative on iiC 1 ), i.e. that we have Vi(x) = i(x), for 
all x G C 1 , where x = dx/dt. As a consequence, the original equation can be recovered 
via the embedded equation by restricting the underlying functional space to i(C k ), k 
depending on the order of the original equation. 

Most of the time, we have in mind applications to physics where Lagrangian systems 
play a fundamental role. The importance of these systems is related to the fact that they 
can be derived via a first principle, the least-action principle. Moreover, in some cases 
of importance we can find a symmetric representation of these systems as Hamiltonian 
systems. Hamiltonian systems are fundamental by many aspects. Of particular interest 
is the fact that they can be quantized in order to obtain quantum analogues of classical 
dynamical systems. 
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An embedding procedure can always be applied to Lagrangian systems. We obtain to 
embed objects which are the embed action functional and the embed Euler-Lagrange equa- 
tion. At this point, an embedding procedure can be considered as a particular quantization 
procedure related to the underlying functional space. However, we have a more elaborate 
picture here. As we have an embed action functional, we can develop the associated cal- 
culus of variations, that we call the embed calculus of variations in what follows. We then 
obtain an embed least-action principle with an associated Euler-Lagrange equation. As a 
consequence, we have two kinds of embed Euler-Lagrange equations and it is not clear that 
the one obtain directly from the embedding of the classical Euler-Lagrange equation and 
the one obtain via the embed least-action principle coincide. This problem is recurrent in 
all embedding theories of Lagrangian systems and is called the coherence problem. This 
problem is far from being trivial in most of the already existing embedding theories, like 
the stochastic one ([18] [20]) or the quantum one [17]. 



In this paper, we develop an embedding theory of ordinary differential equations and 
Lagrangian systems using fractional derivatives. Precise definitions will be given in section 
I. We only give an heuristic introduction for these operators and some basic problems 
where they arise naturally. 

A fractional derivative is an operator which gives a sense to a real power of the classical 
differential operator d/dt, i.e. that we want to consider an expression like 



The previous problem appears for the first time in a letter from Leibniz to L'Hopital in 
1695: " Can the meaning of derivatives with integer order be generalized to derivatives 
with non-integer orders ?" . Many mathematicians have contributed to this topic including 
Leibniz, Liouville, Riemann , etc. We refer to [36] or [32] for a historical survey. A 
number of definitions have emerged over the years including Riemann-Liouville fractional 
derivative, Grunwald-Letnikov fractional derivative, Caputo fractional derivative, etc. 
In this article we restrict our attention to the Riemann-Liouville fractional derivative, 
although the embedding theory can be developed for an arbitrary given fractional calculus 
with different technical difficulties. 



2. Emergence of fractional derivatives 
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The main difficulties when dealing with fractional derivatives are related to the following 
properties: 

(i) fractional differential operators are not local operators 

(ii) the adjoint of a fractional differential operator is not the negative of itself 

Property (i) is widely use in applications and explain part of the interest for these 
operators to model phenomenon with long memory (see for example [12]). 

Other problems arise during computations. Developing the fractional calculus of vari- 
ation and the associated results (the fractional Noether theorem) we have encountered 
difficulties linked with the following facts: 

(i) the classical Leibniz rule (fg)' = f'g + fg' is more complicated (see [40]) 

(ii) there exists no simple formula for the fractional analogue of the chain rule in 
classical differential calculus (see [40]) 

This last difficulty is of special importance in the derivation of the fractional Noether 
theorem. 

The fractional framework has been used in a wide variety of problems. We note in 
particular applications in turbulence [10], chaotic dynamics [46] and quantization [34]. 

In this paper, we will frequently quote the work of F. Riewe ([38], [39]) which proposes 
a fractional approach to nonconservative dynamical systems. The main property of these 
systems is that they induce an arrow of time due to irreversible dissipative effects. The 
relation between fractional derivatives, nonconservative systems and irreversibility have 
been discussed for example in the book [31]. 

Irreversibility implies that we look for the past V t and the future T t of a given given 
dynamical process x(s), s G R at time t, i.e. on the information V t = {x(s),a < s < t} 
and T t = {x(s), t < s < b} where a and b can be chosen and depends on the amount of 
information we are keeping from the past and the future. This induce the fact that we 
look for two quantities, not yet defined that we denote by d-x(t) and d + x(t) from the 
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point of view of derivatives. 

The past and future information can be weighted, i.e. that we look not for x(s) 
but to w(s, t)x(s) where w(s,t) give the importance of the information at time s with 
respect to time t. This can be achieved using a weight uz^s+i an d regularizing the 
corresponding function. We then are lead to two quantities dtx{t) and d+x(t), which 
represent a weighted information on the past and future behavior of the dynamical process. 

The previous idea is well formalized by the left and right (Riemann-Liouville) derivatives 
[40] . We refer to part I for precise definitions. 

3. Deformation theories and the fractional framework 

The fractional framework follows the general strategy outline in the previous section. 
However, a new ingredient comes into play which makes the fractional embedding different 
from the existing stochastic or quantum embedding theories. We have used in this paper 
the left and right Riemann-Liouville derivatives with different indices for the left and right 
differentiation, i.e. we consider a D" and t D^. The extended operator depends naturally 
on these two operators and is denoted V a,/3 . However, this operator does not reduce to 
the ordinary derivative on the set of differentiable functions. We recover the ordinary 
derivative only when a — f3 — 1. As a consequence, we can associate to a given ordinary 
differential equation a two parameters family of fractional differential equations. The 
original equation being recovered for a special choice of these parameters, i.e. a = (3 = 1. 
In that case, we propose to the use the terminology of fractional deformation and to keep 
the terminology of fractional embedding for the procedure which associated a fractional 
analogue of an ordinary differential equation. 

Deformation theory can be formalized as follows: 

A deformation theory is the data of: 

— A finite set V = {(pi, • • ■ ,Pu), Pi £ A} of parameters where A is a given interval of 
R. 

— A //-parameter family of functional spaces T = {J--p}-p & A v ■ 

— Operators T> v defined on T-p such that there exists Vq £ A v satisfying C l C J--p and 
V Vo (x) = x for x £ C 1 . 
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The condition on A is only here to be sure that we have a continuous dependance of the 
whole construction on the parameters. 

The main difference between deformation and embedding lies in the fact that it is no 
usually easy to obtain information on the initial equation from the deformed one. We must 
use asymptotic methods, looking for the behavior of the deformed equation when p — > po. 
This is not the case for a true embedding theory as the initial equation is already present 
in the embedded one. 

4. Plan of the paper 

Our paper has the same architecture as our previous monograph [20] with Sebastien 
Darses about the stochastic embedding of dynamical systems. As a consequence, the 
comparison between the two embedding procedure will be easier. 

In part I we recall the definitions of the left and right fractional derivatives. We also 
define left and right fractional derivatives which have satisfy a semi-group property and 
the adequate functional spaces on which they are defined following a previous work of 
Erwin and Roop [23] . We also recall a product rule formula. 

In part II we define the fractional embedding of differential operators and ordinary 
differential equations. 

In part III we study the fractional embedding of Lagrangian systems. We obtain a 
fractional analogue of the Euler-Lagrange equations. 

In part IV we develop a fractional calculus of variations associated to the fractional 
embedding of classical functionals. generalizing a previous work of O.P Agrawal [1] We 
prove two versions of the least action principle depending on the underlying authorized 
space of variations. We prove in particular a coherence theorem, which roughly speaking 
state that the fractional embedding of the Euler-Lagrange equation coincide with the 
fractional Euler-Lagrange equation obtained via the fractional least-action principle. 

Part V study the behavior of symmetries under the fractional embedding procedure. 
In particular, we prove a fractional Noether theorem which generalizes a recent result of 
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Frederico and Torres [24]. 

In part VI we derive the analogue of the Hamilton formalism for our fractional La- 
grangian systems. 

In part VII we extend results of parts III and IV to cover the Ostrogradski formalism 
for Lagrangian systems. In this unified framework, we recover classical results of F. Riewe 
([38] [39]). Precisely, we obtain a fractional Lagrangian derivation of Nonconservative 
systems. 

In part VIII we study the fractional embedding of continuous Lagrangian systems. 
In particular, we prove that the classical fractional wave equation introduced by W.R. 
Schneider and W. Wyss [41] under ad-hoc assumptions, is the fractional embedding of the 
classical wave equation which respects the underlying continuous Lagrangian structure of 
the equation. An analogous result is obtained for the fractional diffusion equation. 

We then conclude with some open problems and perspectives. 
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PART I 



FRACTIONAL OPERATORS 



In [1] Agrawal has studied Fractional variational problems using the Riemann-Liouville 
derivatives. He notes that even if the initial functional problems only deals with the 
left Riemann-Liouville derivative, the right Riemann-Liouville derivative appears naturally 
during the computations. In this section, we construct an operator combining the left and 
right Riemann-Liouville (RL) derivative. We remind some results concerning functional 
spaces associated to the left and right RL derivative. In particular, we discuss the possibility 
to obtain a law of exponents. 



1.1. Left and Right Riemann-Liouville derivatives. — We define the left and right 
Riemann-Liouville derivatives following ([36] [40] [37] [32]). 

Definition 1.1 (Left Riemann-Liouville Fractional integral) 

Let x be a function defined on (a,b), and a > 0. Then the left Riemann-Liouville 
fractional integral of order a is defined to be 



Definition 1.2 (Right Riemann-Liouville Fractional integral) 

Let x be a function defined on (a,b), and a > 0. Then the right Riemann-Liouville 
fractional integral of order a is defined to be 



Left and right (RL) integrals satisfy some important properties like the semi-group 
property. We refer to [40] for more details. 

Definition 1.3 (Left and Right Riemann-Liouville fractional derivative) 

Let a > 0, the left and right Riemann-Liouville derivative of order a, denoted by a D" 
and t D% respectively, are defined by 



1. Fractional differential operators 



(1.1) 




(1.2) 




(1.3) 



F(n — a) 



1 
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and 

(1.4) tDW0 = _^(_|)"^ (t _ sr «- Ms)ds , 

where n is such that n — 1 < a < n. 

If a = m, m e W, and x G C m (]a, 6[) we have 
(1-5) aD?x=±£, t V? 



dt m ' 1 b dt n 

This last relation which ensures the gluing of the left and right Riemann-Liouville (RL) 
derivative to the classical derivative will be of fundamental importance in what follows. 

d~^~ x d x 

If x(t) G C° with left and right-derivatives at point t denoted by — — and — — respec- 
tively then 



(1-6) a D" 1 x = —, t D™ , . 

v ; 1 dt b dt 

In what follows, we denote by °E, and "E^ the functional spaces defined by 

(1.7) "E={iG C([a, b}), a D?x exists}, Ej = {i G C([a, b}), t D%x exists}, 
and 

(1.8) «Ef = a a En Ef . 

Remark 1.1. - - Of course the set "E^ zs non-empty. Following ([40] Lemma 2.2 p. 35) 
we have AC([a, b]) C "E^, where AC([a,b]) is the set of absolutely continuous functions 
on the interval [a,b] (see [40] Definition 1.2). 

The operators of ordinary differentiation of integer order satisfy a commutativity prop- 
erty and the law of exponents (the semi-group property) i.e. 

d n d m d m d n d n+rn 

(1.9) — o = o — = . 

v ' dt n dt m dt m dt n dt n+m 

These two properties in general fail to be satisfied by the left and right fractional RL 
derivatives. We refer to ([32] §.IV.6) and ([26] p. 233) for more details and examples. These 
bad properties are responsible for several difficulties in the study of fractional differential 
equations. We refer to [37] for more details. 
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1.2. Left and right fractional derivatives. — In some cases, we need that our frac- 
tional operators satisfy additional properties like the semi-group property. Following [23] 
we introduce the left and right fractional derivatives as well as convenient functional spaces 
on which we have the semi-group property. 

Definition 1. 4 (Left fractional derivative). - - Let x be a function defined on R, a > 
0, n be the smallest integer greater than a (n — 1 < a < n), and a = n — a. Then the left 
fractional derivative of order a is defined to be 

jn i mi rt 

(1.10) Wx(t) := ^Dfxit) = ^ooD^xit) = J Jt - s^xtfds. 

Definition 1.5 (Right fractional derivative). - - Let x be a function defined on M., 
a > 0, n be the smallest integer greater than a (n — 1 < a < n), and a = n — a. Then the 
right fractional derivative of order a is defined to be 

(1.11) D?s(f) := t D?x(t) = (-ir— t D^x(t) = jf (* - t^xtfds. 
If Supp(:r) C {a, b) we have Wx = a D?x and B°x = t D%x. 

In [23] several useful functional spaces are introduced. Let / C 1 be an open interval 
(which may be unbounded). We denote by Cq°(I) the set of all functions x G C°°(I) that 
vanish outside a compact subset K of I. 

Definition 1.6 (Left fractional derivative space). Let a > 0. Define the semi- 
norm 

(1-12) | X |ja(R)H| Wx || L 2 (R) , 

and norm 

( \ 1 / 2 

( L1 3) || x ||ja(R):= || x \\ mK) + \x i2 



LV'V Vll uu-yitL) ■ i lj£( 



and let J£(K) denote the closure o/C£°(R) with respect to 



J 



Similarly, we can defined the right fractional derivative space. 

Definition 1.7 (Right fractional derivative space). — Let a > 0. Define the semi- 
norm 

( L14 ) | X |j«(R)H| B * X h 2 (R), 
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and norm 

(2 2 \ 1 ^ 2 

(1-15) || x ||jg(R):= (J| x || i2 ( R ) + I x |ja( M ) J 

and let J#(K) denote the closure of Cq°(R) mt/j respect to \\ • ||jg(R). 

We now assume that / is a bounded open subinterval of R. We restrict the fractional 
derivative spaces to /. 

Definition 1.8. - - Define the spaces J£ (/) ; Jrq{I) as the closure of C£°(I) under their 
respective norms. 

These spaces have very interesting properties with respect to © and D*. In particular, 
we have the following semi- group property: 

Lemma 1.1. — For x G J^ (I), < a < [3 we have 

(1.16) B^x = WB p ~ a x 
and similarly for x G Jr :0 (I), 

(1.17) B%x = B^~ a x. 
We refer to ([23] Lemma 2.9) for a proof. 

The fractional derivative spaces J2,o{I) an d J have been characterized when a > 0. 
We denote by Hq(I) the fractional Sobolev space. 

Theorem 1.1. — Let a > 0. Then the «/£ )0 (/), Jro(I) and H^(I) spaces are equal. 

We refer to ([23] Theorem 2.13) for a proof. In fact, when a 7^ n — 1/2, n G N we have 
a stronger result as the Jl$(I), Jr,o(I) an d Hq(I) spaces have equivalent semi- norms and 



norms. 



2. The extension problem 

As we want to deal with dynamical systems exhibiting the arrow of time, we need to 
consider the operator a D" and t D^, in order to keep track of the past and future of the 
dynamics. The fact that we consider a 7^ (3 is only here for convenience. This can be used 
to take into account a different quantity of information from the past and the future. 

Let a D" and t D^ be given. We look for an operator V a,lS of the form 
(1.18) P a ^ = M( a D t a , t Df), 
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where M : IR 2 — > C is a mapping which does not depends on (a, (3), satisfying the following 
general principles: 

d m x 

— i) Gluing property: If x(t) G C 1 then when a = f3 = m, m G N*, V m ' m x(t) = ——. 

— ii) M is a IR-linear mapping. 

— iii) Reconstruction: The mapping M is invertible. 

Condition i) is fundamental in the embedding framework. It follows that all or con- 
structions can be seen as a continuous two-parameters deformation of the corresponding 
classical one^ . This can be of importance dealing with the fractional quantization problem 
in classical mechanics. 

Condition ii) does not have a particular meaning. This is only the simplest dependence 
of the operator V with respect to a D" and t D^. 

Condition iii) is important. It means that the data of T> a,/3 on a given function x at point 
t allows us to recover the left and right RL derivatives of x at t, so information about x in 
a neighborhood of x(t). 

Lemma 1.2. - - Operators satisfying conditions i), ii) and iii) are of the form 
(1-19) V"* = [p a D? + (p - 1) t Df] + iq [ a Df + t D?], 

where p, q eR and q ^ 0. 

Proof. — By ii), we denote M(x, y) = px + qy + i(rx + sy), with p, q,r,s G K. By i), we 
must have with y = —x corresponding to the operator a choice of operators d/dt, —d/dt 

(1.20) p-q = l, r-s = 0. 

We then already have an operator of the form (1.19). The reconstruction assumption only 
impose that q ^ in (1.19). □ 

A more rigid form for these operators is obtained imposing a symmetry condition. 

( 3 ) Condition i) is not the usual condition underlying the stochastic or quantum embedding theories. In 
general, we have an injective mapping i from the set of diffcrentiable functions C 1 in a bigger functional 
space £ such that the operator V that we define on £ reduces to the classical derivative on /.(C 1 ) meaning 
that for x E C 1 , V{l(x)) = u{x'{t)) where x'(t) = dx/dt. As a consequence, we have a true embedding 
in this case, meaning that the embed theory already contain the classical one via the mapping t. Here, 
the classical theory is not contained in the embedded theory but can be recovered by a continuous two- 
parameters deformation. 
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— iv) Let x G C° be a real valued function possessing left and right classical derivatives 

, , . d + x d_x d + x d_x 

at point t, denoted by — — and —— respectively. It —— = — , then we impose 

dt dt dt dt 

that 

(1.21) V hl x(t) = 

Condition iv) must be seen as the non-differentiable pendant of condition i). Indeed, 
condition i) can be rephrased as follows: if x G C° is such that d + x and d~x exist and 
satisfy d + x = d~x then T> l,1 x = d + x. Condition iv) is then equivalent to the commutativity 
of the following diagram, where R 2 is seen as the R-vector space associate to C: 

C -U c 

x — ix i — > x + ix 

(1.22) MI | M 

C C, 

a i — > ia 

where r : C — > C is defined by r(z) = iz, z G C and we have used the fact that M(x, —x) = 
following condition i). 

Lemma 1.3. - - The unique operator satisfying condition i), ii), Hi) and iv) is given by 
(1-23) V a , p = i^Df - t D?] + <i[ Df + t Df]. 

Proof. — By iv), we must have p + (p — 1) = and 2q — 1, so that p — q — 1/2. □ 



3. The fractional operator of order (a, (3), a > 0, /3 > 

Lemma 1.3 leads us to the following definition of a fractional operator of order (a, (3): 

Definition 1.9. - - For all a,b G R, a < b, the fractional operator of order (a, (3), a > 0, 
(3 > 0, denoted by , is defined by 

(1-24) Vf = \ [ a D? - t Df] + [ a D? + t Df ] , 

where /i G C. 

When a = /3 = 1, we obtain P^'^ = d/di. 

The free parameter /i can be used to reduce the operator T>^ to some special cases of 
importance. Let us denoted by x(t) a given real valued function. 
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- For [i — — i, we have P°'" = a D" then dealing with an operator using the future state 
denoted by J-t{x) of the underlying function, i.e. J-t{x) = {x(s), s e [a, £[}. 

- For jj, — i, we obtain X>°'^ = — then dealing with en operator using the past state 
denoted by Vt(x) of the underlying function, i.e. Vt(x) = {x(s), s &]t,b]}. 

As a consequence, our operator can be used to deal with problems using a D°, t D^, or 
both operators in a unified way, only particularizing the value of /z at the end to recover 
the desired framework. 

When a = — oo and b = oo, we denote the associated operator by D^, i.e. 



as long as f(a) = f(b) = or g(a) = g(b) = 0. 

This formula gives a strong connection between a D" and via a generalized integra- 
tion by part. This relation is responsible for the emergence of i^t m problems of fractional 
calculus of variations only dealing with a D". See section 3 for more details. This result 
also justifies our approach to the construction of a fractional operator which put on the 
same level the left and right RL derivatives. 

As a consequence, we obtain the following formula for our fractional operator: 
Lemma 1.4- — For all f,g E "E^ ; we have 



(1.25) 




where ji G C. 



4. Product rules 



The classical product rule for Riemann-Liouville derivatives is for all a > 



(1.26) 




(1.27) 




provide that f(a) = f(b) = or g(a) = g(b) = 0. 
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Proof. 
(1.28) 



We have 

V a /f(t)g(t)dt 



/(/) 



„ - X ) + MPS 



(g(t))dt. 



Exchanging the role of a and (3 in ( t D° — a Df ) + i//(tD™ + a Df ), we obtain the operator 
( t Df - a D t a ) + i^( t Dj + a D?) which can be written as 

(1-29) - [( a D° - t Dft - i^( t Df + a D°)] = -X>*?. 

This concludes the proof. □ 
Here again, we see that it is convenient to keep the parameter \i free. 
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PART II 

FRACTIONAL EMBEDDING OF DIFFERENTIAL OPERATORS 



1. Fractional embedding of differential operators 

Let d G N be a fixed integer and a, b G M, a < b be given. We denote by C([a, b]) the 
set of continuous functions x : [a, b] — > Let n G N, we denote by C n ([a,b]) the set of 
functions in C([a, 6]) which are differentiable up to order n. 

Let / : M x C d — > C be a function, real valued on real arguments. We denote by F the 
corresponding operator acting on functions x and defined by 

(ni) p . C([a,b}) — > C([a,6]) 

where is the function defined by 

Let f = {fi}i=o,...,n be a finite family of functions, fa : R x C d — > C, and Fi, i = 1, . . . ,n 
the corresponding family of operators. We denote by O f the differential operator defined 
by 

(n-3) °f = E^-^ 

i=0 

where • is the standard product of operators, i.e. if A and B are two operators, we denote 
by A ■ B the operator defined by (A ■ B)(x) = A(x)B(x) and o the usual composition, 

cP 



i.e. (A o B)(x) = A(B(x)), with the convention that y— ) = Id, where Id denotes the 
identity mapping on C. 



Definition II. 1 (Fractional embedding of operators). Let f = {fi}i=o,..., n and 
S = {9i}i=o,...,n be finite families of functions, f i :'RxC d ^C and gi : R x C d — > C 
respectively, and F i} Gi, i = 1, . . . ,n the corresponding families of operators, and Of the 
associated differential operator. 
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The (a, (3) -fractional embedding of Of written as (II. 3), denoted by °Emb^(//)(Of ) is 
defined by 

n 

(II.4) "Emb^XOf) = J> • (VfY o G t . 

i=0 

Note that the embedding procedure acts on operators of a given form and not on 
operators like abstract data, i. e. this is not a mapping on the set of operators. 

We can solve this indeterminacy using a formal representation of an operator. 

Let f = {fi}i=o,...,n and g = {gi}i=o,..., n be finite families of functions, fa : R x C d — > C 
and (?j : K. x C d — > C respectively, and F;, G i: i = 1, . . . ,n the corresponding families of 
operators. We denote by ^Of the operator acting on £ <g> C n <8> £ defined by 

(H.5) Of = ^F,® — 

i=0 

where <g> is the standard tensor product. 

We denote by the set of operators of the form (II. 5) and O the set of differential 
operators of the form (II. 3). We define a mapping n from to O by 

n n 
j=0 i=0 

where /x is the projection n : £ ® £ ® £ — > £, (j,(x <S> y <S> z) — x • (y o z). 

A differential operator being given, its fractional embedding depends on its writing as 
an element of O®. 



2. Fractional embedding of differential equations 

Let k G N be a fixed integer. Let f = {fi}i=o,..., n and g = {gi}i=o,...,n be finite families of 
functions, /,:lx C kd — > C and ft:Rx C kd — > C respectively, and F { , d, i = 1, . . . , n the 
corresponding families of operators. We denote by Of the operator acting on (C n [a, b]) k 
defined by 

n /J* 

(n-7) °? = 5>-^°^ 

i=0 
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The ordinary differential equation associated to Of is defined by 

(11.8) Of(*,^,...,^) = 0, xeC n+k ([a,b]). 

We then define the fractional embedding of equation (II. 8) as follow: 

Definition II. 2. - - The fractional embedding of equation (II. 8) of order (a, (5), a, (3 > 
is defined by 

(11.9) ?Emb?M (Of) (x, Vfx, {Vf) k x) = 0, x G «Ef (n + fc). 

Note that as long as the form of the operator is fixed the fractional embedding procedure 
associates a unique fractional differential equation. 

3. About time-reversible dynamics 

The fractional embedding procedure associate a natural fractional counterpart to a given 
ordinary differential equation. In some case, the underlying equation possesses specific 
properties which have a physical meaning. One of this property is the time-reversible 
character of the dynamics: 

A dynamics on a space U is time-reversible if there exists an invertible map i of U such 
that i 2 = Id, i.e. i is an involution, and if we denote by t the flow describing the dynamics 
we have 

(11.10) io0„ t = <f> t oi, 

meaning that if x(t) is a solution then i(x(—t)) is also a solution of the underlying equation. 

Time- Reversibility is closely related to a specific property of the classical derivative under 
time-reversal: 

(H.H) £ w _ t))= _!<_«), 

We define a notion of reversibility directly on operators: 

Definition II. 3. - - We denote by Rev the C- linear operator defined by Rev( a D") = 
t D?, RevGD") = a D?. 

The action of Rev on is non-trivial: 

fj, 

Lemma II. 1. - - We have Hew(V a /) = -V^. 
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We then have the following analogue of reversibility in the fractional setting: 

Definition II. 4- - - Let Of be a differential operator of the form (II. 7) such that the 
dynamics of the associated differential equation (II. 8) is time-reversible, i.e. 

(11.12) Rev(Of) = -Of. 

The fractional embedding "Embf(/i) is called reversible if 

(11.13) RevCEmbf („) (Of)) = -^Emb^) (Of) . 

The main consequence of lemma II. 1 is that there exists a unique way to do a fractional 
embedding conserving the reversibility symmetry. 

Theorem II. 1. - - The reversibility symmetry is preserved by a fractional embedding if 
and only if a = f3 and /i = 0. 

Proof. — The reversibility symmetry is preserved if and only if we always have 
Rev('D^' /3 ) = — T>® ,/3 . By lemma II. 1 this is only possible when = and a = [3. □ 

In what follows we denote by Rev" the fractional embedding "Emb^(0). 
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PART III 



FRACTIONAL EMBEDDING OF LAGRANGIAN SYSTEMS 



In this section, we derive the fractional embedding of a particular class of ordinary differ- 
ential equations called Euler-Lagrange equations which governs the dynamics of Lagrangian 
systems. 

1. Reminder about Lagrangian systems 

Lagrangian systems play a central role in dynamical systems and physics, in particular 
for classical mechanics. We refer to [3] for more details. 

Definition III.l. — An admissible Lagrangian function L is a function L : 

C such that L(t,x,v) is holomorphic with respect to v, differentiate with respect to x and 

real when oGK. 

A Lagrangian function defines a functional on C 1 (a, b), denoted by 



(III.l) £ a , 6 : C\a,b) -> R, xeC\a,b)< — > L(s, x(s), — (s))ds, a,beR. 



The classical calculus of variations analyzes the behavior of C under small perturbations 
of the initial function x. The main ingredient is a notion of differentiable functional and 
extremals. 

Definition III. 2 (Space of variations). - - We denote by Var(a, b) the set of func- 
tions in C 1 (a, b) such that h(a) = h(b) = 0. 

A functional £ is differentiable at point x G C 1 (a, b) if and only if 



for e > and all h G Var(a, b). 

Using the notion of differentiability for functionals one is lead to consider extremum of 
a given Lagrangian functional. 

Definition III. 3. — An extremal for the functional C is a function x G C 1 (a, b) such 
that dC(x, h) = for all h G Var(a, b). 




(III.2) 



C(x + eh) — L(x) = edC(x, h) + o(e), 
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Extremals of the functional C can be characterized by an ordinary differential equation 
of order 2, called the Euler-Lagrange equation. 

Theorem III.l. - - The extremals of C coincide with the solutions of the Euler-Lagrange 
equation denoted by (EL) and defined by 



(III.3) 



d_ 

di 



dL 

dv 



t,x(t),^(t) 



dL 

dx 



t,x(t),^(t) 



This equation can be seen as the action of the differential operator 

^ dL 



(III.4) 



d dL 

dt dv 



dx 



dx 



on the couple (x(t), ^"(0)- 



2. Fractional Euler-Lagrange equation 

The fractional embedding procedure allows us to define a natural extension of the clas- 
sical Euler-Lagrange equation in the fractional context. The main result of this section 
is: 

Theorem III. 2. — Let L be an admissible Lagrangian function. The "Emb^(/i)- 
fractional Euler-Lagrange equation associated to L is given by 



dL 



dL 



- (t, x(t),V«/x(t)) = — (t, x(t),Vfx(t)) . (FELy 



dx 



In what follows, we will simply speak about the fractional Euler-Lagrange equation 
when there is no confusion on the underlying embedding procedure. 

The proof is based on the following lemma: 

Lemma III.l. - - Let L be an admissible Lagrangian function. The fractional embedding 
of the Euler-Lagrange differential operator O(el) is given by 

dJ dl 

(ni.5) ^M(o iEL) ) = vfo--—. 

Proof. — The operator (III. 4) is first considered as acting on (C 1 ([a, b}) 2 , i.e. for all 
(x(t),y(t)) e C l [a,b] x C^a.b] we have 

(III.6) {EL) (x(t),y(t)) = - ^—( t ,x(t),y(t))j - —(t,x(t),y(t)). 
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This operator is of the form Of with 

dL 



< IIL7 > f= 

and 

(....8) g =(-*l 

where 1 : K. x C 2 — > C is the constant function l(t,x,y) = 1. As a consequence, O^^) is 
given by 

(III.9) o^^.io^-g.Ido!, 



with the convention that = Id. We then obtain equation (III. 5) using definition 

ILL □ 

We can now conclude the proof of theorem III. 2 using definition II. 2. The fractional 
embedding of equation (III. 3) is given by 

(111.10) ?Emb?0i) (0 (EL) ) (x, Vfx) = 0, 

which reduces to equation (FEL^) thanks to lemma III. 1 . 



3. The coherence problem 

The fractional embedding procedure allows us to define a natural fractional analogue of 
the Euler-Lagrange equation. This result is satisfying because the procedure is fixed. How- 
ever, Lagrangian systems possess very special features. In particular, the classical Euler- 
Lagrange equation can be obtained using a variational principle, called the least-action 
principle and denoted LAP. The least action principle asserts that the Euler-Lagrange 
equation characterizes the extremals of a given functional associated to the Lagrangian. 
We then are lead to the following problem: 

— i) Develop a calculus of variation on fractional functionals. 

— ii) State the corresponding fractional least-action principle, in particular explicit the 
associated fractional Euler-Lagrange equation denoted by FEL/; ap . 

— iii) Compare the result with the embedded Euler-Lagrange equation (FEL^ ^) 
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An embedding procedure is called coherent when the two Euler-Lagrange equations are the 
same, i.e. if 

(HI.11) FEL //ap = «Emb£(/i)(£L), 

assuming that FEL// ap is obtained from the embedding of the classical functional using 
the same embedding procedure. 

As we will see, an embedding procedure is not always coherent. Although we obtain in 
general equations of the same form, we usually have some torsion between the embedding 
of the functional and the Euler-Lagrange equation which cancel only in particular cases. 

The fractional calculus of variations is developed in §.IV as well as the corresponding 
fractional least-action principle. The coherence of fractional embedding procedures is dis- 
cussed in §.4. 
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PART IV 

FRACTIONAL CALCULUS OF VARIATIONS 

This section is devoted to the fractional calculus of variations using our fractional op- 
erator. The functional is obtained under the fractional embedding procedure. We refer to 
the work of O.P. Agrawal [1] for related results. 

1. Fractional functional 

Let L be an admissible Lagrangian function on R x R d x C d , d > 1, and C the associated 
functional. Using the Fractional embedding procedure °Embf(/x), we define a natural 
Fractional functional associated to L. 

We denote by "E^ the set of functions x such that a D^x and tD^r are defined. 
Definition IV. 1. - - The Fractional functional associated to L is defined by 
(IV.l) £#:«Ef-R, sG^E?— > [ b L( S ,x( S ),V"/x(s))ds, a,beR; 

J a 

The extension property implies that reduce to the classical functional C a ^ when 
a = P = l. 

2. Space of variations and extremals 

Let us denote by „E^(0, 0) the set of curves h G "E^ satisfying h(a) = h(b) = and 
a E% := a a E%. We denote by Var a (0, 0) the set defined to be 

(IV.2) Var a (a, b) = {h G a K, K a ) = H b ) = and a D?h = t D^h }. 

We denote by P the set 2E£(0, 0) or Var a (a, b). 

Definition IV.2. - - Let x be a given curve. A P '-variation of x is a one-parameter e G R 
family of curves of the form 

(IV.3) y e = x + eh, h G P. 

A notion of differentiability can now be defined for fractional functionals. In the follow- 
ing, we write £ a , b indifferently for C a a f when P = £E£(0, 0) and when P = Var Q (a, b). 
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Definition IV. 3. - - Let L be an admissible Lagrangian function and C a ^ the associated 
fractional functional. The functional C a ^ is called F-differentiable at x if 

(IV.4) £ a , b (x + eh) - C a ,b( x ) = ed £a,b( x i h ) + °( e )> 

for all h G P ; e > ; where dC a ^(x, h) is a linear functional of h. 

The linear functional dC a ,b(x,h) is called the P- differential of the fractional functional 
£ a ,b at point x. 



An extremal for £ a b is then defined by: 

Definition IV. 4- — A F-extremal for the functional C a ^ is a function x such that 
dC a)b {x, h) = for all he P. 

The following lemma gives the explicit form of the differential of a fractional functional: 

Theorem IV. 1. — Let L be an admissible Lagrangian function and the associated 
fractional functional. The functional C"f is differentiable at any x G "E^(x a ,Xb) and for 
all h G ™E^(0, 0) the differential is given by 



(IV.5) d£$(x,h)= f 

J a 



-(t,x(t),V°/x(t)) 



AT 

+ —(t,x(t),V«/x(t)) 



h(t)dt. 



Proof. — As the left and right RL derivatives are linear operators we have 
(IV.6) V% p (x + eh)= Dx + e V a /h. 

As a consequence, we obtain 

(IV.7) C a a f{x + eh)= / L(s,x(s) + eh{s) ,V a / x{s) + eV a /h{s))ds 

J a 

which implies, doing a Taylor expansion of L(s,x(s) + eh(s),V^ ,l3 x(s) + eV^h(s)) in e 

around 

(IV.8) 



— (s,x(s),V»/x(s))h(s) + °—{s,x{s),V«/x{s))Vfh{s) 
+C a a f(x) + o(e). 
Using the product rule (1.27) we obtain 

^(s,x(s),V"/x(s))Vfh(s)ds= - J V^(s,x(s),V a /x(s))h(s)ds. 
Replacing this expression in (IV.8), we deduce formula (IV.5). 



ds 



□ 
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3. The fractional Euler-Lagrange equation 

We obtain the following analogue of the least- action principle in classical Lagrangian 
mechanics: 

Theorem IV. 2 (Fractional least-action principle). Let C[x] be a functional of the 
form 

f b 

(IV.IO) C[x]= / L(s,x(s),V^x(s))ds 

J a 

defined on "E^(£ a ,£&). 

A necessary and sufficient condition for a given function x G "E^ to be a "E^- extremal 
for C[x] with fixed end points x(a) = x a , x(b) = Xb, is that it satisfies the fractional 
Euler-Lagrange equation (FEL): 



(IV. 11) v^l 



dh 

dv 



dL 



(t, x(t),V"/x(t)) = — (t, x(t),V a /x(t)) 



Note that this equation is different from the one obtained via the fractional embedding 
procedure. 

Proof. - - Using the classical Du Bois Reymond lemma ([4], p. 108) and theorem IV. 1 we 
obtain (IV. 11). □ 

The weak analogue using the space of variation Var a (a, b) is given by: 

Theorem IV. 3 (Weak fractional least-action principle). Let C[x] be a func- 
tional of the form 

(IV.12) C[x}= j L(s,x(s),V^x(s))ds 

J a 

defined on a E^{x a ,Xb) . 

A necessary and sufficient condition for a given function x e a E^ to be a Var a (a, b)- 
extremal for C[x] with fixed end points x (a) = x a , x(b) = Xb, is that it satisfies the fractional 
Euler-Lagrange equation FEL^' Q 

We denote FEL£ for FEL£' a in the following. 



4. Coherence 



The coherence problem can now be studied in details. We have the following theorem, 
which is only a rewriting of theorem IV. 2 and definition III. 2: 
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Theorem IV. 4- — Let L be an admissible Lagrangian function, a,b G W, a < b, a, (3 > 0, 
then the following diagram commutes: 

"Embf(u) 

L(t,x(t),dx/dt) L(t,x(t),V^x) 
(IV. 13) LAP | | FLAP 

EL a _&v w FEL7£. 

Theorem IV. 4 is not a coherence result in the spirit of ([15], [16]) or [20]. Indeed, the 
embedding procedure changes between the Euler-Lagrange equation and the functional. 

There exist at least two ways to restore coherence of the embedding procedure. 

The first one is to restrict the set of variations we are looking for in the least-action 
principle. Using the set Var a (a,6) we obtain a coherent procedure via the weak least 
action principle. 

Theorem IV. 5 (Weak coherence). - - Let L be an admissible Lagrangian function, 
a,b G R, a < b, a > 0, then the following diagram commutes: 

Fmb" 

L(t,x(t),dx/dt) a ^ b L(t,x(t),V%x) 
(IV. 14) LAP | | WLAP 

EL FEL" 

where WLAP stands for the weak least-action principle and a Emb^ := "Emb^. 

The second way is to use a specific fractional operator with a particular symmetry 
property. Using the reversible embedding and the operator d a which corresponds to the 
operator T>Q' a , i.e. 

(IV.15) d a := a 1 + t \ 

we obtain the following coherence result: 

Theorem IV. 6 (Reversible coherence). - - Let L be an admissible Lagrangian func- 
tion, a,b G M, a < b, a > 0, then the following diagram commutes: 

L(t,x(t),dx/dt) L{t,x{t),d a x) 
(IV. 16) LAP | I FLAP 

EL ^ FEL a . 

where Rev a is the a-reversible embedding and FEL a denotes the equation FELq'". 
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It is not clear for the moment to know what is the interest of the fractional reversible 
embedding procedure for applications. At least it can be considered as a good candidate 
to a fractional quantization procedure. 
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PART V 

SYMMETRIES AND THE FRACTIONAL NOETHER THEOREM 



In classical mechanics Noether's theorem gives a strong connection between group of 
symmetries of the Lagrangian and conservation laws. We refer in particular to Arnold's 
presentation of Noether's theorem for Lagrangian systems ([3] p. 88). For a historical point 
of view on different generalization of Noether's theorem we refer to [28]. 



1. Invariance of fractional functionals 

Group of symmetries are classical in mechanics. A natural way to deal with symmetries 
is to look for one-parameter group of diffeomorphisms. We first define the action of a 
diffeomorphism on a couple (x(t), V^ ,l3 x(t)). 

Definition V.l. - - Let <p : R d — > R d be a diffeomorphism. The fractional linear tangent 
map associated to (ft, denoted by 0* is defined by 

(V.l) Mx,V a /x) = (<j ) (x),V a /<j ) (x)). 

Definition V.2. - - Let $ = {0 s } se R be a one-parameter family of diffeomorphisms. An 
admissible Lagrangian L is said to be invariant under the action of $ if 

(V.2) Hx,V a /x) = L(0 s (x),P^(0 s (z))), G R, 

or equivalently 

(V.3) L(x, V a /x) = L ((<j>,Ux, V a /x)) . 

The relation with infinitesimal transformations used in [24], is obtained using a Taylor 
expansion of yt(s) = <f> s (x(t)) in a neighborhood of 0, leading to 

(V.4) y t (s)=y t (0) + s.^(0) + o(s). 

As O = Id is the identity map, we obtain denoting -^(0) = £(t,x) an infinitesimal 

as 

transformation given by 

(V.5) x(t)t-+x(t) + s£(t,x(t)) + o(s). 

A natural problem in the embedding framework is the following: assume that L is 
invariant under a one-parameter group of diffeomorphisms. What can we said about the 
fractional embedded Lagrangian ? 
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2. The fractional Noether theorem 



In this paragraph, we derive a fractional version of the Noether theorem. Similar results 
have been obtained by G. Frederico and D. Torres [24] in a different setting. 

Theorem V.l (The fractional Noether theorem). Let L be an admissible La- 
grangian invariant under a one-parameter group of diffeomorphism $ = {</> s } se R. Then, 
we have 



(V.6) 



o ( <9L\ dy t dh p f dy t \ _ 
\0v J ~ds ~dv {its') |s= °-°' 



where yt(s) = (f> s (x(t)), along all solutions of the fractional Euler-Lagrange equation. 

Equation (V.6) is not the usual way to state the Noether theorem in classical mechanics. 
This is mainly due to the fact that the fractional differential operator T>^ ,lS inherits from the 
bad properties of the underlying left and right Riemann-Liouville derivatives concerning 
the Leibniz rule. 

Proof. — Let yt(s) = (p s (x(t)), the £"b invariance is equivalent to 



(V.7) 



-[L(t,y t (s),V«/y t (s))]=0. 



The usual chain rule for the classical derivative implies 

(V.8) ^(t,yt(s),V a /y t (s)) ■ ^ + ^(t,y t (s),V a /y t (s)) ■ ± [Vfy t (s)} = 0. 



ds dv 



ds 

dy t (s) 



As P"" 9 acts on the variable t and d/ds on the variable s, and — ; — G "Ef, we deduce 
M ds 



that 
(V.9) 

As a consequence, we obtain 
dL 



d 



dyt 
ds 



dx (t,y t (s),V"/y t (s)) ■ ^f- + ^(t,y t (s),V«/y t (s)) • Vf 



(V.10) 

As x(t) is an extremal for L, we have 
(V.ll) 

This concludes the proof. 



dyt 
ds 



0. 



\dv~ 



dL 

dx 



□ 
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3. Toward fractional integrability and conservation laws 

Of course, when a — (5 — 1, equation (V.6) reduces to 
( , d (dh\dy t dhd (dy t \ _ 

(V ' 12) It lb + TvJt \-d~s) Is= °- °' 

which is equivalent to 

<™> a (§>)-* 

using the Leibniz rule. 

In that case, the Noether theorem can then be stated in term of conservation laws (first 
integrals), i.e. saying that the quantity 

(V.14) C(t,x(t),x(t)) = ^(t,x(t),x(t)) ■ ^(0), 

is constant along all the solutions of the Euler-Lagrange equation. 

First integrals play a fundamental role in classical Lagrangian dynamics as they are 
related to the classical problem of integrability (Liouville's integrability), i.e. the fact 
that a Lagrangian systems with sufficiently many first integrals can be integrated by 
quadratures (see [3] Chapter 10). 



It is not clear regarding to the fractional Noether's theorem to define the analogue of 
conservation laws and in particular the fractional analogue of integrability for Lagrangian 
systems. 



JACKY CRESSON 



FRACTIONAL EMBEDDING 



37 



PART VI 

FRACTIONAL HAMILTONIAN SYSTEMS 



The classical passage from the configuration space, i.e. space of positions to the phase 
space, i.e. the space of position and velocities, is achieved using the Hamiltonian formalism. 
In this section, we derive the fractional analogue of Hamiltonian systems. 



1. Hamiltonian systems and Legendre property 

Let L be a Lagrangian system. When the Lagrangian satisfies an analytic property, 
that we called the Legendre property, we can find a remarkable symmetric formulation of 
the Euler-Lagrange equation called Hamiltonian. The Legendre property can be stated as 
follows: 

Definition VI. 1. - - Let L be an admissible Lagrangian function. We say that L satisfies 

dL 

the Legendre property if the mapping v i— > — (x, v) is invertible for all x G M. d . 

As a consequence, there exists a map / : M. d x C d — > C d such that for any p G C d 
dL 

satisfying p = —(x,v) we have v = f(x,p). The map / is called the Legendre transform. 

The Legendre transform connects the momentum 
(VI.l) p(t) = —(x(t),x(t)), 

of a given solution x(t) of the Euler-Lagrange equation with its velocity x(t). 

Using the Legendre transform, we can define a new function called the Hamiltonian 
associated to L. 

Definition VI.2. - - Let L be an admissible Lagrangian system satisfying the Legendre 
property, and f the corresponding Legendre transform. The Hamiltonian function associ- 
ated to L, denoted by H : M. d x C d — > C 7 is defined by 

(VI.2) H(x,p)=pf(x,p)-L(x,f(x,p)). 
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The Hamiltonian system associated to H is given as 

dx dH 

< VL3 > I I %H 

dt dx ' 

The dynamics of the Hamiltonian system is equivalent to the dynamics governed by the 
Lagrangian. 



2. The fractional momentum and Hamiltonian 

The previous construction can be carried in the fractional context for our fractional 
Euler-Lagrange equation. 

Definition VI. 3. - - The fractional momentum of a given solution x(t) of the fractional 
Euler-Lagrange equation is defined as 

(VI.4) p{t) = °—{x{t)^/x(t)). 

We remark that the fractional momentum is obtain from the classical one by the frac- 
tional embedding procedure. If L satisfies the Legendre property we obtain a fractional 
analogue of the relation between momentum and velocities: 

Lemma VI. 1. - - Let L be an admissible Lagrangian system satisfying the Legendre prop- 
erty and f the associated Legendre transform. Let x(t) be a solution of the fractional 
Euler-Lagrange equation, and p(t) its fractional momentum. We have 

(VI.5) V a /x(t) = f(x(t),p(t)). 

The fractional Euler-Lagrange equation is then equivalent to the following system of 
fractional differential equations: 

V a /x(t) = f(x(t),p(t)), 

(VL6) V^p{t) = ^(x(t),f(x(t),p(t))), 

which is nothing else that 

f)T~f 

V«/x{t) = —f(x(t),p(t)), 
(VL7) % 

Vifrit) = - — (x(t),p(t)), 
with H the Hamiltonian function (VI. 2). 



We then introduce the following definition of a fractional Hamiltonian system: 
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Definition VI. 4- - - Let L be an admissible Lagrangian system satisfying the Legendre 
property. The fractional Hamiltonian system associated to L is defined by equation (VI. 7). 



3. Fractional Hamilton least-action principle 

Fractional Hamilton equations (VI. 7) can also be obtained via a variational principle 
called fractional Hamilton leat-action principle. 

Theorem VI. 1. - - We denote by H : M. d x C d — > C a Hamiltonian function and F the 
associated functional defined by 

(VI.8) F(x(t), p(t)) = f {p{t)Vfx{t) - H(x(t), p(f))) dt. 

J a 

A couple (x(t),p(t)) such that x(a) = x a , p(a) = p a , x(b) = Xb, p(b) = pb is an extremal 
of (VI.8) if and only if it satisfies the fractional Hamiltonian equations (VI. 7). 

Proof. - - We denote by L the function 

(VI. 9) L(x, p, v, w) — pv — H(x, p). 

The functional F can be seen as a functional associated to the Lagrangian L. The fractional 
Euler-Lagrange equation associated to this functional is 

(Vi-io) g W()) ) = 

(vi.il) ■*$ gw*))) - 

(VI. 12) 

where z(t) = (x(t),p(t),Vfx(t),Vfp(t)). 

d\j ,9L dH , n . 

As — — = p and — — = — — — , the first equation gives 
ov ox ox 

PjTT 

(VI.13) V^p(t) = - — (x(t),p(t)). 

Moreover, as 7— = and — — — v — — , the second equation gives 
ow op op 

(VI.14) = V"/x(t)-—(x(t),p(t))- 

This concludes the proof. □ 
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As a consequence, the usual construction of Hamiltonian equations using the Legendre 
transform and fractional momentum is coherent with the fractional Hamilton least-action 
principle. 
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PART VII 



FRACTIONAL OSTROGRADSKI FORMALISM AND 
NONCONSERVATIVE DYNAMICAL SYSTEMS 



This part is devoted to the fractional embedding of functionals depending on higher 
order derivatives, i. e. to a fractional analogue of the Ostrogradski formalism. We derive 
the generalized fractional Euler-Lagrange equation. As an application, we prove that some 
Nonconservative dynamical systems can be obtained via the fractional embedding of well 
chosen higher order Lagrangian systems. In particular, we obtain a new approach to results 
obtain previously by F. Riewe ([38] [39]). 



We first introduce some convenient definitions and notations. 

Definition VII. 1. - - Let n G N* and d G N* ; a Lagrangian function of order n denoted 
by L is a function defined by 



forx G C d , Vi G C d , i = l,...,n. 

A Lagrangian function of order n possesses a natural Lagrangian functional denoted by 
C and defined by 



1. Ostrogradski formalism and fractional embedding 



(VII.l) 



R x C d x C d x • • • x C d 
. ' 

n times 

(t,X,Vx, ...,v n ) 



+ L(t,x,ui,...,u n ), 



C, 




J a 



where x^ n \t) denotes the n-th derivative of x with respect to t, i.e. 



(VII.3) 




where 



(VII.4) 




n times 



with o the usual composition of operators. 



42 



JACKY CRESSON 



The associated notion of extremals differs from definition III. 3 only on the boundary 
conditions, i. e. on the space of variations. 

Definition VII. 2 (Space of variations of order n). — Let n G W , we denote by 
Var(a, b) n the subset of functions h G C n (a,b) satisfying 

h(a) = h'(a) = ... = h^- 1 \a) = 0, 
{V 0) h(b) = h'(b) = ■■■ = h^(b) = 0. 

We deduce the following generalized notion of extremals. 

Definition VII. 3. - - Let n G N* , L be a Lagrangian function of order n and C the 
associated functional. The functional C is called differentiate at x G C n (a,b) if and only 
if 

(VII.6) C(x + eh) - C(x) = edC(x, h) + o(e), 

for all h G Var(a, b) n where dC(x, h) is a linear functional of h. 

We then define extremals of the functional C as follows: 

Definition VII. 4- — An extremal for the functional £ is a function x G C n (a,b) such 
that dC(x, h) = for all x G Var(a, b) n . 

Extremals of the functional C are characterized by a generalized Euler-Lagrange equation 
given by 

'dL 



^(tMtU 1 \t),...^(t)) + J2(-iy 

i=l ^ ' 



^(t,x(t),xU(t),...,x^(t)) 

OVi 



(GEL) 

We refer to ([25] Chap. 2 §.11 p. 42) for a proof. We denote by OLAP the Ostrogradski 
least-action principle, i.e. the derivation of (GEL) from C. 

The generalized Euler-Lagrange equation can be written using a differential operator 
denoted by Ogel and defined by 

The (GEL) is then equivalent to 

(VII.8) GEL .(t,:r(t),:r«(t), . . .,x<- n \t)) = 0. 
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The fractional embedding procedure can be applied to obtain a fractional analogue of 
this generalized Euler-Lagrange equation. As for the fractional Euler-Lagrange equation, 
we need to impose some regularity properties on the Lagrangian. 

Definition VII. 5. - - Let nGf and L(t, x,v±, . . . , v n ) be a Lagrangian function of order 
n. The Lagrangian L is admissible ifL is differentiate with respect to x, holomorphic in 
each variable v i} i — 1, . . . , n, and its restriction tolxl^x (R d ) n is real. 

We have: 

Theorem VII. 1. - - Let n G N* and L be an admissible Lagrangian function of order n. 
The "Emb^(/x) -fractional embedding of the generalized Euler-Lagrange equation associated 
to L is given by 



BT n 

—(t,x(t),x^(t),...,x^(t))+ j £(-iy wy 

i=l 



^(t,x(t),x^(t),...,x^(t)) 



a,0 



where 

(VII.9) x°®f> = (T% >f, y x, i = l,...,n. 

The proof is based on the following lemma: 
Lemma VII. 1. - - Let n e N* and L be an admissible Lagrangian function of order n. 



■is 



The "Emb^(/x) -fractional embedding of the generalized Euler-Lagrange operator Ogel 
given by 

(VII.10) -EmbfM (Ogel) = B" 1 )* W*j ° ^ + W 

i=i * 

Proof. — The operator (VII. 8) is of the form Of with: 

- f = {fi}i=Q,..., n where = (-1)* and (— l)*(t, x,v 1 ,..., v n ) = (-1)* for % = 0, . . . , n. 

- g = {9i}i=o,...,n where g = ^ and g { = i = l,...,n. 

n ■ d i 

We then have G el = ^^tF l ■ -z— o t Gi- Using the fractional embedding "Embf(/x) of 

(XL 

i=0 

operators (II. 4) we obtain (VII.10). □ 

By definition II. 2, the fractional embedding of equation (VII. 8) is given by 
(VII.11) ^Embf (jm) (Ogel) • (s, E^s, . . . , (P^) n ) = 0, 
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which reduces to (FGEL) thanks to lemma VII. 1. This concludes the proof of theorem 
VII. 1. 

As for the fractional Euler-Lagrange equation, we have an associated coherence problem 
(see §.3). Next section deals with the corresponding fractional calculus of variations. 



2. Fractional Ostrogradski formalism 

We develop the notion of extremals and differentiability for functionals depending on 
higher order, i.e. power of the fractional operator . We prove in particular a coherence 
theorem. 

2.1. Space of variations and extremals. — The fractional analogue of the space of 
variations of order n is: 

Definition VII. 6 (Fractional space of variations of order n) 

Let n EW , we denote by "Var^(n) the set of functions h E "E^ such that 

(VII. 12) (V^yh{a) = (V^Yh(b) =0 for % = 0, . . . , n - 1. 

2.2. The fractional generalized Euler-Lagrange equation. — The explicit form of 
the differential of a fractional functional of order n E W is given by the following lemma: 

Theorem VII. 2. - - Let JiGf and L be an admissible Lagrangian function of order n 
with £ its associated functional. The functional £ is differentiable at any x E "E^(n) and 
for h E "Varf (n) the differential is given by 



(VII.13) 



d£(x, h) 



V 



dL 

dvi 



z n {t)) 



dL 

dx 



(*»(*)) 



dt, 



where z n (t) = (t, x(t),V a /x(t), (V a /) n x(t)) . 

Proof. — By linearity of we have V^(x + eh) = V^x + eV^h. We deduce that 

,6 

(VII.14) £(x + eh) = / L(z n (s) + e(0,h(s),V^h(s),...,(V^) n h(s)))ds. 

J a 

Doing a Taylor expansion of L(z n (s) + e(0, h(s),V^ ,l3 h(s), . . . , (V^) n h(s))) in e around 0, 
we obtain 
(VII. 15) 

£(x + eh) = 



dL 



dL 



Hzn(s)) + e[J2 [(Pf) ih ^)\ ^«( s )) + Ks) — {z n {s)) ) + 0(6 



,i=i 



dx 



ds. 
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Using the product rule (1.27), we obtain by induction the formula 

(vine) f f(t)(v«/yh(t)dt = [\-iy \(v^yf(t)] h(t)dt, i > 1. 

Replacing in equation (VII. 15) we obtain (VII. 13). □ 

We obtain the fractional analogue of the least-action principle in the fractional Ostro- 
gradski formalism, which follows easily from theorem VII. 2. 

Theorem VII. 3 (Fractional Ostrogradski least-action principle) 

Let L be an admissible Lagrangian of order n G N* 7 and C its associated functional. A 
necessary and sufficient condition for a given function x G "E^(n) to be an extremal for £ 
with fixed end points is that it satisfies the fractional generalized Euler-Lagrange equation 
(FGEL)t$. 

We denote by FOLAP the previous derivation of (FGEL). 

2.3. Coherence. — As for the fractional embedding of classical mechanics, the fractional 
embedding of the Ostrogradski formalism is not always coherent. We have the usual torsion 
on the exponents and the parameter //. 

Theorem VII. 4- - - Let L be an admissible Lagrangian of order n, a, b G R, a, (3 > 0, 
then the following diagram commutes 

T f . d/dtf / 1 \ \ \ a Emb fe(/^) T f . V«'P f M u 
L(t,Z n ' {X{t))) > L(t,Z n ^ {X(t))) 

(VII. 17) OLAP| | FOLAP 

GEL fEmb^) pGEL ^ 

where z®(x(t)) = (x(t), 0(x(t)), . . . , O n (x(t)), with O an operator which can be d/dt or 

Here again, in order to obtain a coherent diagram, we must choose /i — and a — (5 at 
the end, i.e. on the diagram (VII. 17). 



3. Nonconservative dynamical systems 

In this section, we provide a unified framework for F. Riewe's approach to noncon- 
servative dynamical systems ([38], [39]) using the fractional embedding procedure. For a 
different approach to this problem we refer to the seminal paper of P.J. Morrison [33]. 
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3.1. Linear friction. — We consider the classical problem of linear friction 

,„ rTT _ „ s dU dx d 2 x 

< VIU8) -* + ^ + m *7 = °- 

where 7 G i, m > 0. 

This system can not be modelled by an ordinary Lagrangian systems. This follows from 
Bauer theorem [9], which states that it is impossible to use a variational principle to derive 
a single linear dissipative equation of motion with constant coefficients. 

We begin with the following immediate lemma. 
Theorem VII. 5. - - Let L:KxCxCxC^Cfcfte Lagrangian of order 2 defined by 
(VII. 19) L(t,x, Vl ,v 2 ) = -U(x) - ^vf + l -mv\, 

where 7 6 I, m > and U a smooth function. 

The D°- fractional embedding of the generalized Euler- Lagrange equation associated to L 
is given by 

(VII.20) -^ + 7 (D^) 2 :r + m(D^) 4 a; = 0. 

Specializing to /i — — i, we obtain D" i = D a . The operator D a satisfies a semi-group 
property (see lemma 1.1). As a consequence, we obtain WB? = D a+/3 for a > 0, (3 > 0. 
We deduce the following corollary: 

Corollary VII. 1. - - The D° -fractional embedding of the generalized Euler- Lagrange 
equation associated to L is given by 

(VII.21) -^ + 7 B 2a i + mD 4a i = 0. 

dx 

For a = 1/2 we obtain only fractional operators of integer order, which reduce to ordinary 
derivatives on the set of smooth functions. As a consequence, we have proved: 

Theorem VII. 6. - - The D 1 / 2 -fractional embedding of the generalized Euler- Lagrange 

equation associated to L reduce to the equation 

, T TTT nnS dU dx d 2 x 

on the set of smooth functions. 

Using the weak-coherence lemma we deduce the stronger result: 
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Theorem VII. 7. - - The solution of the non-conservative systems (VII. 18) corresponds 
to smooth weak-extremals of the © 1//2 -fractional functional associated to L. 

This lemma is not valid for a general Lagrangian systems. However, in our particular 
case, we have a very simple dependance of the Lagrangian with respect to velocities. This 
is this property with the fact that the restriction on C 2 -functions kills the imaginary part, 
which implies the result. 

3.2. The Whittaker equation. — It is in general difficult to see directly if an equation 
is Lagrangian or not. Following Bateman [5], E.T. Whittaker proposes as an answer to a 
problem raised by R.C. Tolman ([43], §.10 §.11) that the equation 

(VII.23) I -H I 

called the Whittaker equation in what follows, can not be derived by any Lagrangian 
approach. It seems that up to the author knowledge the question is still open today. 

Following a previous work of F. Riewe [39] we give a positive answer to R.C. Tolman's 
problem in the fractional framework. 

We begin with the following immediate lemma. 

Lemma VII. 2. - - Let L : R x C 6 be a Lagrangian of order 2 on C 2 defined by 

(VII.24) L(t, x, y, vi, ui, v 2 , u 2 ) =v\ + u\- u 1 v 1 + x 2 + v 2 y. 

The D°- embedding of the Euler-Lagrange equation associated to L is given by 

rV TT9^ 2x(t)+D«(D^)-(D«) 2 (2(D°)Mt)+y(t)) = 0, 

1 j i(0+B«(D«x)-(D«) 2 (2(D«) 2 i/(0) = 0- 

Specializing when \x = —i, we obtain D™ i = 3 a . The operator D Q satisfies a semi-group 
property (see lemma 1.1). As a consequence, we have D Q D a = B> 2a for all a > and 
equation (VII. 25) reduces to: 

Theorem VII. 8. - - The B) a -embedding of the Euler-Lagrange equation associated to L is 
given by 

2x(t)+3 2a y-2B 4a x(t)-I} 2a y(t) = 0, 
lViL2bj x(t) + B) 2a x - 2B ia y(t) = 0. 
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Choosing a = 1/2, we obtain fractional operators of integer order. On the set of suffi- 
ciently smooth functions, these operators reduces to the ordinary derivatives. As a conse- 
quence, we have proved: 

Theorem VII. 9. - - The restriction of the O 1 / 2 -fractional embedding of the Euler- 
Lagrange equation associated to L is given by the Whittaker equation. 

Moreover, using the weak least action principle, we have the stronger result: 

Theorem VII. 10. — Solutions of the Whittaker equation correspond to smooth weak- 
extremals of the 1//2 -fractional functional associated to L. 
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PART VIII 

FRATIONAL EMBEDDING OF CONTINUOUS LAGRANGIAN 

SYSTEMS 

The wave equation is an example of a partial differential equation which can be derived 
via a variational principle. In that case, the functional depends on several variables (at least 
two). Such functional arise in mechanical problems involving systems with infinitely many 
degrees of freedom like strings, membrane, etc. These systems are known as continuous 
Lagrangian systems. We use the fractional embedding procedure to define the fractional 
analogue of these equations. In particular, we obtain a fractional Lagrangian formulation 
of the fractional wave equation introduced by Schneider and Wyss [41]. This suggest also a 
canonical way to generalize to the fractional framework a given differential equation: if the 
underlying equation possesses an additional structure like the Lagrangian case, one must 
provide an extension which respect this structure. The fractional embedding procedure 
is well adapted to deal with such an extension. We also derive a fractional Lagrangian 
formulation for the fractional diffusion equation. 

1. Continuous Lagrangian systems 

Let d e N. We consider a Lagrangian function L defined on R x R d x R x C x R d and 
denoted by 

(VIII.l) li(t,x u . . . . . .,w d ). 

In what follows, we use the terminology of Lagrangian density for a function L of the form 
(VIII.l). We denote x for (x±, . . . 

Let 9ft be a fixed region of R d and a < b, a,b e R. We consider the functional 



(VIII.2) 




acting on a function u : 1 x t 
in all its variables and where 



R which is usually called a field, which is of class C 1 
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with h satisfying the boundaries conditions 

(VIII.5) h(a, x) = h(b, x) = and h(t, 93?) = 0, 

where <93? denotes the boundary of 3?. 

Theorem VIII. 1 (Euler- Lagrange equation for fields). Extremals of the func- 
tional (VIII. 2) are solutions of the Euler- Lagrange equation for fields (ELF): 



(V1II.6) 



dL d_ fdL 

dy dt \dv 



i=i 



_d_ 

dxi 



dL 



0. 



We refer to ([25] Chapter 7) for more details. 



As in the previous parts, we associate to (ELF) the differential operator 

?n f) at . d 

(VIII. 7) 0/ 



d dL 

o 



(ELF) Qy q v Z_/ g x . q w . ' 

acting on 2d + 2-uplet (t, x, u(t, x), d t u(t, x), d Xl u(t, x), . . . , d Xd u(t, x)). 



Using the fractional embedding procedure we obtain a fractional Euler-Lagrange equa- 
tion for fields. 

Theorem VIII. 2. — Let L be an admissible Lagrangian density. The "Embf (/i)- 
fractional embedding of the Euler-Lagrange equation for fields (ELF) associated to L is 
given by 



AT 

-{z a , p {t,x))-Vf 



dL, . ' 

-g^{z at p{t,x) 



d 

y— 



1=1 



dL 

dwi 



(FELF)Z 



a,/3 



where z a> p(t,x) = (t, x,u(t, x),V^u(t, x), d xi u(t, x), . . . , d Xd u(t, x)). 



Proof. — The differential operator O(elf) is considered as a time-differential operator, i.e. 
that we consider the operator Of with 



(VIII.8) 
and 

(VIII.9) 



— - V — o — 1 

dy ^ dxi dwi ' 



g 



dL\ 



Using definition II. 2, we obtain (FEL). This concludes the proof. 



□ 
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The previous theorem can be generalized in various ways. In particular, one can define 
a generalized embedding procedure, assuming that not only time derivatives, but also 
spatial derivatives d/dxi, % — 1, . . . , d, are replaced by a fractional derivatives. In this case 
however, we have some technical difficulties (see remark VIII. 1). This will be discussed in 
another work. 



As for the previous cases, we develop the associated fractional calculus of variations. 



2. Fractional continuous lagrangian systems 

Using the fractional embedding procedure, we look for the following class of fractional 
densities: 

Definition VIII. 1. - - Let L be an admissible Lagrangien density. The fractional func- 
tional associated to L is defined by 

(VIII.10) C a a f{u)= f I ' L(t,x,u(t,x),V^u(t,x),d x u(t,x)) dxdt, 

for fields u(t, x) G "F^(K) ; the set of fields smooth with respect to x and in "E^ with respect 
to t. 

We consider two spaces of variations for fields: 

Definition VIII. 2 (Spaces of variations for fields). We denote by Var a /3 5R (a, 6) 

the set of fields satisfying 

(VIII.11) 

Var Q ' /3 (a, b, M) = {h(t, x),h t eC 1 , h x G ?e£, h(a, x) = h(b, x) = 0, h(t, dU) = 0}, 
and by Varg(a,6) the set of fields defined by 

rVTTT 11) Vxv a (n h sftt - / M^)A <E C\ h x G JEf, h(a,x) = h(b,x) =0,1 

As usual the condition a D"^ = tD& h allows us to obtain a symmetric product rule which 
is fundamental for the coherence problem. 

Definition VIII. 3. - - Letu(t,x) be a given field. A ^-variation of u is a one-parameter 
e G R family of fields u t (t,x) defined by 

(VIII.13) u e (t,x)=u(t,x) + eh(t,x), he¥, 

where P can be Var a ' /3 (a,6) orVar£(a,6). 
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We denote by C afi the functional C a a f when P = Var a>/3 (a, b, 9fc) and when P = 
Var Q (a,6, 9fc). 

Definition VIII. 4- - - Let L &e an admissible Lagrangian density and C a ,b the associated 
functional. The functional C a ^ is called P- differentiate at u, where u is a field, if 

(VIII.14) C a)b (u + eh) - C a>b (u) = ed£ a , b {u, h) + o(e), 

/or a// h G P, where dC a ^(u, h) is a linear functional of h. 

The differential of the fractional functional (VIII. 10) is given by: 

Theorem VIII. 3. - - Let L be an admissible Lagrangian density and C^f the associated 
fractional functional. The functional is differentiate at any fields uF and for all 
h G Var Q,/3 (a, 6, 3?) the differential is given by 



(VIII.15) QjjMt,x))-V^ 



dh, . ' 



Ed 
ftnr- 



i=l 



dL 

dwi 



z a ,p{t,x)) 



where z a>/3 (t,x) = (t, x,u(t, x),V^u(t, x), d xi u(t, x), . . . ,d Xd u(t,x)). 

Proof. — As the left and right (RL) derivatives are linear operators, we have 

(VIII. 16) V a /(u + eh) = V a /u + eV a /h. 

As a consequence, we obtain 

(VIII.17) C°f(u + eh)= I I L(z u (t,x) + ez h (t,x)) dxdt, 

Ja JU 

where for a field f(t,x) we denote 

(VIII.18) z f (t, x) = (t, x, f(t, x),V a /f(t, x), d xi f(t, x),..., d Xd f(t, x)). 

Doing a Taylor expansion of L(z u (t,x) + ezh(t,x)) in e around 0, we obtain 
(VIII. 19) 



L(z u (t,x) + ez h (t,x)) 



L(z u (t,x)) 



dy 
+o(e). 

Using the product rule (1.27) we obtain 

V a /h—{z u {t )X ))dxdt = 



dL dL dL 

— (z u {t, x))z h (t) + V a /h—(z u (t, x)) + -^M^ x )) d *i h 



duii 



(VIII.20) 



dv 



hV 



dL . 



dxdt 



as long as h(a, x) = h(b, x) = 0. 
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Moreover, using the multidimensional Green theorem (see [44] p. 223), we have 

b » d rvr „?, „ d 



as long as h(t, 93?) = 0. This concludes the proof. 



dwi 



z u (t,x)) 



h dx dt, 



□ 



Remark VIII. 1. — In the global fractional framework, i.e. when we use fractional 
derivatives also in the x variables, we see from formula (VIII. 21) that we will need an 
analogue of the multidimensional Green theorem in the fractional case. 

We denote by D; the operator V«' a and F a (3?) the set £F£(&). 

Theorem VIII. 4 (Real variations). - - Let L be an admissible Lagrangian density and 
£-a'b the associated fractional functional. The functional is differentiable at any fields 
u G F a (3?) and for all h G Var a (a, b, 3?) the differential is given by 



(VIII.22) 



dy 



(z a (t,x))-V" 



dh 

dv 



(z a (t,x) 



i=i 



_d_ 

dx,. 



dL 



z a (t,x)) 



where z a (t, x) = (t, x, u(t, x),V^u(t, x), d Xl u(t, x), . . . , d Xd u(t, x)). 



The proof follows the same lines. 



An immediate consequence of this lemma is: 



Theorem VIII. 5 (Fractional least-action principle for fields) 

Let L be an admissible Lagrangian density and the associated functional. A neces- 
sary and sufficient condition for a field u to be a Var Q,/3 (a, b, ^t)-extremal is that it satisfies 
the fractional Euler- Lagrange equation for fields (FELF)^ . 

We denote by the functional C^t- Using the set of real variations, we have: 

Theorem VIII. 6 (Weak Fractional least-action principle for fields) 

Let L be an admissible Lagrangian density and C^ b the associated functional. A necessary 
and sufficient condition for a field u to be a Var Q (a, 65ft) - extremal is that it satisfies the 
fractional Euler-Lagrange equation for fields (FELF)^. 
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3. The fractional wave equation 

In this section, we derive the fractional wave equation defined by Schneider and Wyss 
[41] as the extremals of a fractional continuous Lagrangian systems. 

The equation describing waves propagating on a stretched string of constant linear mass 
density p under constant tension T is 

rVTTT OK d 2 u(t,x) (PU^X) 

(VIIL23) P ^r~ = ' 

where u{t,x) denotes the amplitude of the wave at position x along the string at time t. 
The wave equation corresponds to the extremals of the generalized functional associated 
to the Lagrangian systems 

(VIII.24) L(t,x,^, W ) = |t; 2 -- W 2 . 

In [41], the authors define the fractional analogue of the wave equation by changing the 
classical derivative by a fractional one. Using our notations, the definition of the fractional 
wave equation is: 

Definition VIII. 5. - - The fractional wave equation of order a is the fractional differ- 
ential equation 

(VIII.25) —pB) 2a u = T-. 

A natural demand with respect to this generalization which is just a formal manipulation 
on equations, is to keep a more structural property of the wave equation, namely the fact 
that it derives from a least-action principle. Using our fractional embedding procedure, we 
are able to explicit such a fractional Lagrangian framework for the fractional wave equation. 

In the following we work with the fractional embedding associated to which corre- 
sponds to V"' 13 when (5 — a and a = — oo, b = +oo. 

Theorem VIII. 7. - - The 3^-fractional embedding of the continuous Euler- Lagrange 
equation associated to (VIII.24) is given by 

(VIII.26) -pD« o D> = Td 2 x2 u. 

We can specialized by choosing p, — —%. In that case = D Q and satisfies a semi-group 
property (see lemma 1.1). As a consequence, we obtain: 
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Corollary VIII. 1. - - The U) a -fractional embedding of the continuous Euler- Lagrange 
equation associated to (VIII. 24) is given by 



Moreover, using the weak coherence theorem, we have: 

Theorem VIII. 8. - - Solutions of the fractional wave equation (VIII. 27) of order a cor- 
responds to weak-extremals of the D Q -fractional functional associated to L. 

Up to the author knowledge, this is the first time that the fractional wave equation is 
derived via a fractional variational principle. In particular, the previous derivation has 
the advantage to keep the continuous Lagrangian structure underlying the classical wave 
equation. 



The previous remark can be uses as a conceptual guideline to generalize classical equa- 
tions of physics in the fractional framework. If the classical equation possesses an addi- 
tional structure, for example Lagrangian, then we must extend this equation keeping this 
additional structure, generalized in a natural way. The main remark is that equations by 
themselves do not have a universal significance their form depending mostly on the coordi- 
nates systems being used to derive them. On the contrary the underlying first principle like 
the least-action principle carry an information which is of physical interest and not related 
to the coordinates system which is used. At least this point of view explain the importance 
of coherence theorems in all the existing embedding theories of dynamical systems. 



(VIII.27) 



pB 2a u = Td 2 x2 u. 



4. A remark on fractional/classical equations 



5. The fractional diffusion equation 




The aim of this section is to derive a fractional Lagrangian formulation of the fractional 
wave equation of order < a < 1, then including the classical diffusion equation. In the 
contrary of the fractional wave equation, the diffusion equation is recovered thanks to a 
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still fractional variational principle. 

Let us consider the Lagrangian function L defined on R x R 2 x C x C by 

1 a 2 

(VIII.29) L(t, x, y, v, w) = -v 2 - -w 2 , 

where p e R, a G R. 

As for the fractional wave equation, we denote by the quantity when (5 = a, 
a = — oo and b = +oo, i.e. working with the fractional derivatives D anf O*. 

Theorem VIII. 9. - - The D^ 2 -fractional embedding of the continuous Euler- Lagrange 
equation associated to (VIII.29) is given by 

,d 2 u 



(VIII.30) B a J 2 oB a J 2 u = a 2 - ,. 

Choosing ji = —i, we obtain = D a . As D a satisfies a semi-group property, we finally 
obtain: 

Theorem VIII. 10. - - The O"/ 2 '-fractional embedding of the continuous Euler- Lagrange 
equation associated to (VIII.29) is given by 

3 2 u 

(VIII.31) 3 a u = a 2 — . 

ox 2 

It must be noted that even for a = 1, the diffusion equation is recovered using a 
fractional embedding procedure, namely the D 1//2 -fractional embedding procedure. 

The main result of this section is that this fractional embedding of the diffusion equation 
has an additional structure, a Lagrangian one. 

Theorem VIII. 11. - - Solutions of the fractional wave equation (VIII. 28) of order 
< a < 1 corresponds to weak-extremals of the IB)"/ 2 -fractional functional associated to 
(VIII.29). 

This result seems new, even for the case a = 1. 
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CONCLUSION AND PERSPECTIVES 

Part of this paper can be generalized in various ways. 

i) We must extend the fractional embedding procedure in order to cover partial differ- 
ential equations. In particular, the general fractional reaction-diffusion equation must be 
studied in this setting. 

ii) Of special interest is the Schrodinger equation. With F. Ben Adda ([7], [8]) we 
have obtain this equation following an idea of L. Nottale [35] and results about the local 
fractional calculus [6]. A problem is to try to do the same computations using left and 
right RL fractional derivatives. The main difficulty is the absence of a chain rule for these 
operators. In that case, one must obtain a one-parameter family of partial differential 
equations, depending on the non-integer order a of differentiation. Following [16] the idea 
is to prove that the Schrodinger equation can be recovered only for a = 1/2 which gives 
strong constraints on the underlying nature of space-time in Nottale's framework of the 
Scale relativity. 

iii) The previous idea can in fact be generalized to stochastic processes following our 
previous work with S. Darses ([18] [20]). In that case, one want to consider fractional 
Brownian processes via the stochastic embedding procedure. The stochastic derivative 
defined in [20] must be generalized and the paper by S. Darses and I. Nourdin [22] can be 
considered as a first step in this direction. 

iv) We have use the left and right RL fractional derivatives, but the same scheme can 
be used with another fractional calculus leading to a different fractional embedding theory. 
This depends mostly on the type of applications one want to consider. Of particular interest 
seems to be the Caputo fractional derivatives [40]. A forthcoming paper will be devoted 
to this case. 
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NOTATIONS 



Operators 

— a D": left Riemann-Liouville fractional derivative of order a 

— t ^>b : right Riemann-Liouville fractional derivative of order a 

— 3 a : left fractional derivative of order a 

— right fractional derivative of order a 

— T>^\ fractional operator of order (a, (5) and parameter jj, 

— D^ 3 : fractional operator of order (a, (3) and parameter [i associated to D and ©* 

— d a : the reversible fractional derivative corresponding to 2)"^ for (5 = a, /i = 0. 

Functional spaces 

— "E: left fractional Riemann-Liouville derivative space 

— E^: right fractional Riemann-Liouville derivative space 

— "E^: functional space associated to the fractional operator D"" 9 

— AC([a,b\): absolutely continuous functions on [a,b] 

— C^°(I): functions of C°°{I) vanishing ouside a compact subset K of /. 

— J£(R): left fractional derivative space 

— J^(1R): right fractional derivative space 

— J2,o(I), Jr,o(I) : closure of C™(I) under their respective norms 

— Hq(I): fractional Sobolev space of order a 

Operators and objects associated to the fractional embedding procedure 

— "Emb^(/i): the fractional embedding procedure associated to T>^ 

— Rev": the reversible fractional embedding procedure, i.e. "Emb^(0) 
~~ 0(el) : the Euler-Lagrange operator 

— EL: the Euler-Lagrange equation 

— FEL^: the fractional Euler-Lagrange equation 

— Var a : the real space of variations 

— Ca'b'- the fractional functional associated to L 
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